Elimination of transverse instability in stripe solitons by one-dimensional lattices 
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We demonstrate theoretically and experimentally that the transverse instability of coherent soliton 
stripes can be greatly suppressed or totally eliminated when the soliton stripes propagate in a one- 
dimensional photonic lattice under self-defocusing nonlinearity. 



It is well known that in homogeneous nonlinear opti- 
cal media, a bright soliton stripe, uniform along the 
transverse stripe (say y) direction but localized along 
the orthogonal transverse (say x) direction, is unstable 
upon propagation along the longitudinal z direction when 
transverse perturbations are present [l|-l9j. When a one- 
dimensional (ID) optical lattice is introduced along the 
x or y direction, the soliton stripe is still transversely un- 
stable under self- focusing nonlinearity pH [H| • To sup- 
press this transverse instability, some ideas have been 
proposed. For instance, this instability can be completely 
eliminated if the soliton stripe is made sufficiently inco- 
herent along the transverse direction [12j. This instabil- 
ity can also be significantly reduced by nonlinearity sat- 
uration or incoherent mode coupling [l3|, [l4| . The intro- 
duction of 2D square lattices can also suppress transverse 
instability, but the transversely stable structures in such 
media are soliton trains (comprising an infinite array of 
intensity peaks) rather than soliton stripes [HI, [HI . 

In this paper, we demonstrate, both theoretically and 
experimentally, that transverse instability of coherent 
soliton stripes is greatly suppressed or totally eliminated 
when the soliton stripes propagate in a ID lattice under 
self-defocusing nonlinearity. 

Our theoretical model is the 2D NLS equation with a 
ID lattice, 



iU z + U x 



yy 



n(x)U + a\U\ 2 U = 0, (1) 



where a = ±1 denotes self- focusing and self-defocusing 
nonlinearity, and n(x) is a ID lattice. For definiteness, 
we take 

n(x) — —6 sin 2 x 

in this paper (see Fig. Ufa)). Stripe (ID) solitons in this 
model are of the form U (x, z) = u(x)e~ l ^ z , where u(x) 
is a real- valued localized function, and /i is the propaga- 
tion constant. When a = — 1 (defocusing nonlinearity), 
there is a stripe-soliton family in the first bandgap of the 
lattice. The power curve of this family is shown in Fig. 
□Jb). Here, the power P is defined as \u\ 2 dx. The 
intensity profile of the soliton at fi = 4.5 is displayed in 
Fig. [TJc) (the peak intensity is approximately 3.2). To 
determine the transverse stability of these stripe solitons, 



we perturb them as 

U(x,y,z) = e~^ z [u{x) + [v(x) + w(x)]t 



,iky+\z 



+ [v*(x) -w*(x)]e- iky + x * z } 



where w,u]<l are normal-mode perturbations, and k is 
the transverse wavenumber. Substituting this perturbed 
solution into (pQ) and neglecting higher order terms of 
(y,w), we obtain the linear-stability eigenvalue problem 



LqW 



where 



d x 



■ n(x) 

■ n(x) 



k 2 



i\w, 



3au 2 , 



and A is the eigenvalue. The full spectrum of this eigen- 
value problem can be obtained numerically by the Fourier 
collocation method [9]. For the stripe soliton at \i = 4.5 
(see Fig. [U^c)), its stability spectrum is displayed in 
Fig. (He). This spectrum lies entirely on the imaginary 
axis, indicating that this stripe soliton is transversely 
stable under defocusing nonlinearity! In contrast, when 
the nonlinearity is self- focusing (a = 1), stripe solitons 
will remain transversely unstable, and this instability is 
strong. To demonstrate, a family of stripe solitons in the 
semi-infinite gap under focusing nonlinearity is obtained, 
and its power curve is plotted in Fig. [ljb). At \i = 1 
of the power curve, the soliton is shown in Fig. DJd) 
(its peak intensity is roughly 1.7). The linear-stability 
spectrum of this soliton is displayed in Fig. DJf). This 
spectrum contains large positive eigenvalues (with max- 
imum 1.25), indicating strong instability. 

Next we consider the linear stability of other stripe 
solitons in the solution families of Fig. QJb). When the 
solitons are near an edge /io of a Bloch band, these soli- 
tons under perturbations are low-amplitude Bloch-wave 
packets, 

U(x, y, z) = e~^ z [etf (X, Y, Z)p(x) + e 2 U 2 + • • • ] , 

where p(x) is the Bloch wave at edge /io, < e < 1, 
X = ex, Y = ey, Z = e 2 z, and ^(X, Y, Z) satisfies 



i^z + D^xx + ^yy + era |*| 2 * =0, 



(2) 
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FIG. 1: (a) ID lattice; (b) power curves of stripe solitons 
(dashed red indicates instability, solid blue indicates stability, 
and shaded regions are Bloch bands); (c) intensity profile of 
a stripe soliton in the first gap (at ji = 4.5) under defocusing 
nonlinearity; (d) intensity profile of a stripe soliton in the 
semi-infinite gap (at fi = 1) under focusing nonlinearity; (e,f) 
stability spectra of stripe solitons in (c,d) respectively. 



disappear, hence stripe solitons in this /i range are all 
stable! To demonstrate, the most unstable eigenvalue 
Amax versus \i is plotted in Fig. [2j It is seen that the 
real part of Amax (the maximum growth rate) is under 
0.03 for all stripe solitons, indicating that the transverse 
instability is extremely weak. In addition, when \i > 4.3, 
the real part of Amax becomes zero, hence those stripe 
solitons are linearly stable. These results show that under 
defocusing nonlinearity, transverse instability of stripe 
solitons is either greatly suppressed or totally eliminated. 
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FIG. 2: The most unstable eigenvalue Amax versus \i for 
stripe solitons in the first gap under defocusing nonlinearity: 
(left) real part; (right) imaginary part. 



with D being the diffraction coefficient of the ID lattice 
at edge /io and a > being a constant Stripe soli- 
tons in Eq. (pp) correspond to stripe envelope solutions 
\f(X,y,Z) = A(X)e- irZ , where sgn(r) = -sgn(a) = 
— sgn(D), and A(X) is a sech function [9[. It is well 
known that this stripe envelope solution is transversely 
unstable in the envelope equation (|2]) [3,0, HI- Thus low- 
amplitude soliton stripes in Eq. (pQ) are all transversely 
unstable. In particular, for the stripe-soliton family in 
the semi-infinite gap, D > 0, thus the transverse insta- 
bility is of neck- type (due to positive eigenvalues); while 
for the soliton family in the first gap, D < 0, thus the 
transverse instability is of snake-type (due to both pos- 
itive and complex eigenvalues) [9[. The magnitude of 
these unstable eigenvalues is proportional to — /io| 

Away from band edges, we have tracked these 
transverse-instability eigenvalues for the two soliton fam- 
ilies in Fig. QJb). We find that for the soliton family in 
the semi-infinite gap (under focusing nonlinearity), as \i 
moves away (decreases) from the band edge /io = 2.06, 
the magnitude of the largest positive eigenvalue keeps 
increasing (at \i = 1, this magnitude has reached 1.25, 
see Fig. DJf)). Thus all stripe solitons in this family are 
transversely unstable. This instability is of neck- type and 
is strong when \i is not near the band edge. 

However, for the soliton family in the first gap under 
defocusing nonlinearity, the story is very different. In this 
case, as /i moves away (increases) from the band edge 
/io = 2.27, the positive eigenvalues of snake instability 
quickly disappear when \i > 2.3. Meanwhile, the real 
parts of the complex eigenvalues in the snake instability 
first saturate quickly to a very low level when 2.4 < /i < 
4.3. After \i > 4.3, these complex eigenvalues then totally 



The above linear-stability results are also corroborated 
by nonlinear-evolution simulations of these stripe soli- 
tons under random-noise perturbations. For the soliton 
in Fig. []Jc), its initial perturbed state (with 2% random- 
noise perturbations) is shown in Fig. [3fa), and evolution 
output of this perturbed soliton under defocusing non- 
linearity at z = 100 is shown in Fig. [3fb). It is seen 
that even after such a long evolution, this stripe soliton 
still remains robust and does not break up. In contrast, 
when the soliton in Fig. DJd) is perturbed by the same 
amount of perturbations, after evolution under focusing 
nonlinearity, this stripe quickly breaks up into filaments 
at z = 3 (see Fig. [3^c)). Thus suppression of transverse 
instability under defocusing nonlinearity and persistence 
of transverse instability under focusing nonlinearity hold 
for nonlinear evolutions as well. 
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FIG. 3: (a) Initial intensity pattern of the stripe soliton in 
Fig. [TJc) under 2% perturbations; (b) output intensity of 
the perturbed soliton in (a) after nonlinear evolution of z — 
100 under defocusing nonlinearity; (c) output intensity of the 
perturbed soliton in Fig. W[d) after nonlinear evolution of 
z = 3 under focusing nonlinearity. 

Experimentally, we have confirmed the above theoret- 
ical predictions. The experiments were performed in a 
lOmm-long biased SBN:60 photorefractive crystal. The 
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FIG. 4: Experimental results. Upper first column: ID lat- 
tice; lower first column: initial probe beam; upper row: re- 
sults with lattice; lower row: results without lattice; second 
column: linear diffraction; third column: output with focus- 
ing nonlinearity; fourth column: output with defocusing non- 
linearity. The bias field and crystalline c-axis are along the 
vertical direction. 



optically induced ID lattice (41/im lattice spacing) is 
shown in the upper first panel of Fig. [4j while the ini- 
tial stripe beam (12/im FWHM) is shown in the lower 
first panel. The peak intensity ratio between the probe 
and lattice beams is about 1:5. In the presence of the 
lattice, this probe beam exhibits strong discrete diffrac- 
tion during linear propagation (upper second panel). It 
breaks up due to strong transverse instability under self- 
focusing nonlinearity at a bias field of 2kV/cm (upper 
third panel), but remains robust against transverse per- 
turbations under defocusing nonlinearity at a bias field 
of — 1.6kV/cm (upper fourth panel). As a comparison, 
when no lattice is induced inside the crystal, the corre- 
sponding results are shown in the lower panels of Fig. 
[4j In this case, the stripe beam breaks up strongly due 
to neck-type transverse instability under focusing nonlin- 
earity Q , and broadens even more than linear diffraction 
under defocusing nonlinearity. 

It is important to notice from Fig. QJb) and Fig. [2] 
that transversely stable stripe solitons (under defocusing 
nonlinearity) are located near the second Bloch band, 
implying that their stability is caused by mode coupling 
between the first and second bands. This means the ex- 



istence of these stable stripe solitons cannot be predicted 
by the corresponding discrete NLS model, 

lU n ,z + Un+1 ~ 2Z7 n + C/n-i + U n , yy + a\U n \ 2 U n = 0, 

since this discrete model is derived under a single-band 
approximation and it does not incorporate mode cou- 
pling between different Bloch bands [17]. To confirm 
this, we seek stripe solitons in this discrete model as 
U n (y, z) = u n e~ l ^ z , where \i is the propagation constant. 
The power curve of these solitons P = ^2\u n \ 2 versus \i 
for a = — 1 (defocusing nonlinearity) is plotted in Fig. 
[5fa), and the most unstable linear-stability eigenvalue of 
each soliton versus \i is plotted in Fig. [5fb,c). It is seen 
that when \i moves away (increases) from the band edge 
/io = 4, the real part of the most unstable eigenvalue 
quickly saturates (to about 0.52) and then stays at this 
level for all higher /i, thus transverse instability persists 
for all these discrete stripe solitons under defocusing non- 
linearity. 
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FIG. 5: (a) Power curve of discrete stripe solitons under de- 
focusing nonlinearity; (b,c) real and imaginary parts of the 
most unstable eigenvalue Amax versus \i. 

In summary, we have demonstrated both theoretically 
and experimentally that the transverse instability of co- 
herent soliton stripes is totally eliminated or greatly 
suppressed when the soliton stripes propagate in a ID 
photonic lattice under self-defocusing nonlinearity. This 
elimination of transverse instability makes stripe solitons 
applicable in physical settings. 

This work was supported by the NSF and AFOSR. 



[1] V.E. Zakharov and A.M. Rubenchik, Sov. Phys. JETP 
38, 494 (1974). 

[2] N.N. Akhmediev, V.I. Korneev, and R.F. Nabiev, Opt. 

Lett. 17, 393-395 (1992). 
[3] A.V. Mamaev, M. Saffman and A. A. Zozulya, Europhys. 

Lett. 35, 25 (1996). 
[4] Z. Chen, M. Segev, T.H. Coskun, and D.N. 

Christodoulides, Opt. Lett. 21, 1436 (1996). 
[5] A. De Rossi, S. Trillo, A.V. Buryak, and Y.S. Kivshar, 

Phys. Rev. E 56, 4959 (1997). 
[6] A. De Rossi, S. Trillo, A.V. Buryak, and Y.S. Kivshar, 

Opt. Lett. 22, 868-870 (1997). 
[7] Y.S. Kivshar and D.E. Pelinovsky, Phys. Rep. 331, 117 

(2000). 

[8] S.P. Gorza, Ph. Emplit, and M. Haelterman, Opt. Lett. 



31, 1280 (2006). 
[9] J. Yang, Nonlinear Waves in Integrable and Noninte- 

grable Systems (SIAM, Philadelphia, 2010). 
[10] A.B. Aceves, C. De Angelis, G.G. Luther, A.M. 

Rubenchik, Opt. Lett. 19, 1186 (1994). 
[11] D. Neshev, A. A. Sukhorukov, Y.S. Kivshar, and W. Kro- 

likowski, Opt. Lett. 29, 259-261 (2004). 
[12] C. Anastassiou, M. Soljacic, M. Segev, E.D. Eugenieva, 

D.N. Christodoulides, D. Kip, Z.H. Musslimani, and J. P. 

Torres, Phys. Rev. Lett. 85, 4888 (2000). 
[13] Z.H. Musslimani, M. Segev, A. Nepomnyashchy, and Y.S. 

Kivshar, Phys. Rev. E 60, R1170 (1999). 
[14] Z.H. Musslimani and J. Yang, Opt. Lett. 26, 1981 (2001). 
[15] X. Wang, Z. Chen, J. Wang and J. Yang, Phys. Rev. 

Lett. 99, 243901 (2007). 



4 



[16] J. Yang, Phys. Rev. A. 84, 033840 (2011). Equation (Springer, Berlin, 2009). 

[17] P.G. Kevrekidis, The Discrete Nonlinear Schrddinger 



